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Conformal Mapping on Rough Boundaries 
I: Applications to harmonic problems 
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The aim of this study is to analyze the properties of harmonic fields in the vicinity of rough 
boundaries where either a constant potential or a zero flux is imposed, while a constant field is 
prescribed at an infinite distance from this boundary. We introduce a conformal mapping technique 
that is tailored to this problem in two dimensions. An efficient algorithm is introduced to compute 
the conformal map for arbitrarily chosen boundaries. Harmonic fields can then simply be read 
from the conformal map. We discuss applications to 'equivalent' smooth interfaces. We study the 
correlations between the topography and the field at the surface. Finally we apply the conformal 
map to the computation of inhomogeneous harmonic fields such as the derivation of Green function 
for localized fiux on the surface of a rough boundary. 

PACS numbers: 2.70.-C, 66.10.Cb, 44.30.+V, 61.43.Hv 
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Defining and computing effective properties of hetero- 
geneous media is a subject which has been studied for a 
long time, and for which a number of powerful techniques 
have been developped. In most cases however, the hetero- 
geneities are considered to lie in the bulk of the material. 
Another type of inhomogeneity is due to the random ge- 
ometry of the surface on which boundary conditions are 
applied. This study focusses on this second type. We 
will thus consider homogeneous media which are limited 
by a rough surface or interface. Our purpose here is to 
introduce a very efficient way of solving harmonic prob- 
lems in two-dimensional systems for any geometry of the 
boundary. 

The occurence of rough interfaces in nature is more the 
general rule than the exception. Apart from very specific 
cases such as mica where a careful cleavage can produce 
planar surfaces at the atomic scale, surfaces are rough. 
Even glass with a very homogeneous composition, where 
the surface is obtained by a slow cooling of the material, 
so that surface tension can act effectively to smoothen all 
irregularities, does display roughness in the range 5 to 50 
nanometers over a window of a few micrometers width 1]. 
Similarly, the so-called "mirror" fracture surface which is 
optically smooth exhibits specific topographic patterns 
when examined with an atomic force microscope|3. The 
key question is thus to identify the relevant range of scales 
at which roughness appears. From common observations, 
this question may not admit a simple clear cut answer. 
Indeed in a variety of cases, the amplitude of the rough- 
ness appears to be strongly dependent on the size of the 
examined surface. A particular class of such scale depen- 
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dent roughness, namely self-affine roughness^\^, has re- 
cently motivated a lot of activity (see References |M IE Q 
for recent reviews) both because of its relevance in many 
different instances, and of their theoretical justification 
which has been obtained in statistical physics for a wide 
class of models ranging from growth models^], molecular 
beam epitaxy ,§1, fracture surfacesIlOj, to immiscible fluid 
interfaces 11]. Although the present study is not specific 
to self-affinc surfaces, we shall consider this particular 
class in order to apply our method. The interest of this 
choice being that i) the description of the roughness is 
realistic for a number of applications, ii) consequences 
can be expressed in quite general terms as a function of 
few parameters directly accessible experimentally, and fi- 
nally Hi) the most commonly studied roughness models 
are "monochromatic" surfaces with a single asperity pat- 
tern repeated periodically, and hence the transposition 
to more complex geometries may reveal wrong (exam- 
ples of such cases will be explicited in the main body of 
this article). 

As previously mentioned, if most surfaces are rough, 
this roughness may be of small amplitude macroscopi- 
cally, and thus one may feel that its role can be neglected 
in most cases. Fortunately, this is generally true. Tak- 
ing into account precisely the surface roughness may be 
required in two distinct classes of problems: 

The first class (I) covers applications where the rough- 
ness cannot be neglected at the scale at which the bulk 
field varies. For obvious reasons, there is no way to avoid 
the accurate description of the boundary. We may men- 
tion the following potential applications: 
• In confined geometries, such as naturally encountered 
in surface force study, the roughness of the surface may 
affect the interpretation and thus the precision of the 
measurements since the distance between two facing sur- 
faces is generally estimated from indirect measurements 
of transport in the gap between the surfaces, 
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• Field which are rapidly varying in space will be sensi- 
tive to fine details of the boundary geometry. The most 
obvious example in this field is the reflection and scat- 
tering of a wave by a rough boundary. 'l3] Of particular 
importance is the case of surface waves, evanescent waves, 
Rayleigh waves in elasticity, ... 

• In a similar spirit, diffusion processes may display 
anomalous behaviors at short times where the diffusion 
length is smaller or comparable to the roughness. [T^ 

The second class of problems (II) where roughness can- 
not be neglected is when one has to focus on the bound- 
ary, either because only this part matters for extraneous 
reasons or because the system is sensitive to high fields 
which can be induced by the roughness itself. Some ex- 
amples of these two cases are listed below: 

• Surface phenomenon such as electrofiltration requires 
a proper solution of say a Stokes flow field, in the imme- 
diate vicinity (typically Debye length scale) of a rough 
boundary where an electric boundary layer is present and 
can be entrained by the fluid to give rise to an electric cur- 
rent in response to a fluid flow in a porous medium, fill 

• The brittle fracture of glass is generally due to surface 
defects which induce locally high stresses, which reduce 
significantly the breaking limit of this material. In the 
absence of specific surface degradation the most impor- 
tant source of surface defect is the topography itself 16]. 

• Some growth models have a local growth rate which 
depends on a harmonic field locally. The development of 
unstable modes which will finally induce a macroscopic 
roughening does require the proper analysis of the field 
at the surface. 0, 111 

The relative independence of the bulk field on the small 
scale roughness of the boundary for a slowly varying field 
(class II problems) can be used to explore the local field 
close to the boundary using an asymptotic analysis with 
a double scale technique. The large scale problem con- 
sists in solving the problem at hand replacing the rough 
boundary by a smooth equivalent one. The small scale 
problem deals with the details of the rough boundary 
and matches at "infinity" with a homogeneous field. This 
local problem will be considered in full details in the fol- 
lowing. 

These examples are obviously not exhaustive. Inho- 
mogeneous boundary conditions may arise for instance 
in contact problems where the roughness cannot be 
neglected. 17,18] One may also consider application out- 
side the realm of physical applications, such as the use of 
harmonic problems and particularly conformal map for 
providing a simple means of meshing a domain limited 
by a rough boundary. 

In the present article we will essentially focus on har- 
monic problems. The latter arise in a variety of differ- 
ent domains in physics, such as electrostatics, thermal 
or concentration diffusion, flow in porous media, anti- 
plane elasticity to mention a few. Another use of con- 
formal mappings is the resolution of bi-harmonic prob- 
lems near a rough interface; both stress field in elas- 
ticity and velocity field in low Reynolds number fluid 



mechanics IH, IM Ell can be derived from po- 
tentials that obey such bi-laplacian equations. We refer 
the reader to the companion paper which is completely 
devoted to this specific problem. 

This paper is devoted to the study of harmonic prob- 
lems in 2D semi-infinite media limited by a rough bound- 
ary. To extend the definition of the profile of the bound- 
ary to infinity, we use periodic boundary conditions along 
the boundary. Although very specific, this type of geom- 
etry will be very convenient as soon as no other boundary 
lies close to the first one. The distance threshold to con- 
sider in such a case is typically of the order of magnitude 
of the larger spatial wavelength of the profile i.e. the spa- 
tial period in the geometry we have described. We use a 
conformal mapping technique. It consists in contructing 
a map from the domain of interest (in the complex plane) 
onto a regular semi-plane. The conformity of the map al- 
lows to preserve harmonicity through the map transform. 

In a first part of this paper, we define the form of the 
conformal mapping suited to our geometry. Then, we 
adress the question of constructing the mapping associ- 
ated with any prescribed interface. We show that this 
problem can be solved with an iterative algorithm using 
Fast Fourier Transforms (FFT). This algorithm allows 
to get the conformal map in a few iterations of FFT, 
whose computation time scales as Nlog{N) where N is 
the number of Fourier modes used to describe the in- 
terface. Note the remarquable efficiency of such a tech- 
nique, considering that the map gives the solution of a 
Laplacian field in the entire two-dimensional problem. 
This problem is very close to the so-called "Theodorsen 
problem" [2^ in a circular geometry. We also show that 
one can generate maps which naturally give rise to self- 
affine boundaries, a powerful technique to explore generic 
properties of such problems. Specific applications of this 
technique to self-affine profile are studied, which includes 
i) the question of defining an equivalent smooth (planar) 
interface, and finding its height compared to the geomet- 
rical average height of the interface, ii) the correlation 
between the height and the field which is computed ex- 
actly in the limit of a small roughness amplitude. These 
two examples demonstrate the unexpected difference in 
behaviors for persistent and anti-persistent profiles. Fi- 
nally, we give the expression of the Green function for 
localized flux on a rough interface. 



I. TRANSFORMATION CONFORME 

A. Notations 

In order to study harmonic fields in two dimensions, 
very powerful techniques have been developed based on 
complex analysis _29j. Among these, we will use in the fol- 
lowing conformal maps, which allow to relate the geome- 
try we wish to study (i.e. a semi-infinite domain limited 
by a rough interface), to a regular one as schematically 
illustrated in Figure (Q. 
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FIG. 1: Illustration de la transformation conforme Q intro- 
duite dans cette section. A schematic illustration of the map- 
ping fl which maps the semi-infinite plane T> onto the domain 
limited by a rough interface £. 



As usual, we will identify a point in the plane {x,y), 
with the complex number z — x + iy. We note z — x — iy 
the complex conjugate of z. The two variables z and 
z can be treated as independent variables instead of x 
and y. A mapping from the complex plane onto itself 
is simply defined as a complex function O of z and z 
which transforms one point of the complex plane z into 
another point Q{z,z). For the mapping to be of physical 
interest it has to be bijective in a domain of interest, and 
thus inversible. The mapping is conformal if the function 
is holomorphic, i.e. it only depends on z and not on 
J. It can be shown that in this case, local angles are 
preserved in the transformation — apart from singular 
points — and hence the term "conformal" . Moreover, the 
real part and the imaginary part of any such holomorphic 
function are both harmonic, i.e. SRe $7 = 3m H, — 
0. The latter property results from the expression of the 
Laplacian operator in terms of the variables z and z: 



yy zz 



(1) 



Two obvious choices are well suited for our geometry: 
i) a semi-infinite plane V, 5m(z) < ii) the unit circle 
C, \z\ < 1. These two domains can be related by the 
transformation z log(z), and thus they are basically 
equivalent. Since the boundary we consider is periodic 
in the a;-direction, the mapping to the unit circle is well 
suited. However, in the following, we will rather use the 
mapping to the half-plane, since it corresponds directly 
to the "reference" problem where the roughness vanishes. 



B. Specific transformations 

The domain of interest, denoted by £ is limited by a 
rough interface d£ which is a periodic function of 5Re(z) 
of period X. The conformal map, Q, is a function of 
z which associate one point of the reference domain, V 
or C, to another point in £. From now on in order to 
distinguish the initial and the image domain, we will note 
w = u + iv a point in the image plane, and keep the 
notation z = x + iy for the initial plane unless otherwise 
mentioned. Before specifying the particular form of the 
boundary, it is possible to guess an adequate form for 
these transformations. 



Let us first consider the mapping from the half-plane 
V to £. As Sm(z) tends to — cx), the mapping should 
approach the identity, z, since the roughness of 

the boundary is not expected to play any significant role 
at a large distance (compared to the period X) from the 
boundary. We introduce the function uj{z) such that 



n{z) = z + uj{z) . 



(2) 



Functions of the form exp(— i/cz) with k real, thus ap- 
pears to be natural candidates for lo{z). They are indeed 
periodic functions of 3fie(z), and vanish exponentially as 
3m(z) goes to — oo when fc > 0. Moreover, in order 
to satisfy the same periodicity as d£ we require that 
kX ~ 2m: where n is an integer. Thus we propose as 
a writing of the transformation 17 the following decom- 
position 



il{z) = Z + WfcC 
r!=0 



-2i7rfcz/X 



(3) 



Without loss of generality we will set X — 2tt for all the 
rest of this study. 

The rough boundary is to be identified with the image 
of the X axis, so that d£ obeys the parametric equations 



u = X- 
V — 



The corresponding transformation from the unit disk C 
to the domain £ can be obtained from the above form ^ 
and the tranformation from the disk to the semi-infinite 
plane 'ismz < 0. The resulting transformation reads 



il(z) = i log(z) + ^ LUkz'' 
k 



(5) 



where X = 1 has been used. Evidently, the image of the 
unit circle z = e*^ provides the same parametric form as 

O- 

The form of the transformation being imposed, one 
needs to check that the transformation is bijective: a 
point should have a single image, and an image point a 
unique parent. This condition does impose some restric- 
tion on the transformation Q. It can be rephrased simply 
for the transformation Q as 

dn 



dz 



> , (6) 

z=x+iy 

for all y < 0. In principle, it is sufficient to impose this 
condition only in the strict interior of the domain. If 
dVl/dz = on the boundary, a kink may appear at this 
point. In the following, we will assume that the interface 
is smooth at a small scale so that poles are forbidden on 
the boundary. 



As an example, if the transformation is simply 
Jl(z) = z -t- we"*^ 



(7) 



then the condition JB)) reduces to |1 -I- iwe"e*'| > or 
\uj\ < 1. For this maximum value, the image of the x- 
axis is a cycloid, with cusp points. Figure (|2Jl illustrates 
this limit case. 
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FIG. 2: Image of a regular square grid of the semi-plane "D 
using the transformation Q. The parameter u has been set 
to its maximum value of 1. Cusp appears on the boundary. 



II. COMPUTING THE MAPPING FOR AN 
IMPOSED INTERFACE 

The above presented transformation is only useful for 
a particular application if the transformation f2 can be 
computed, once the boundary d£ is imposed. This sec- 
tion is devoted to this problem. The algorithm that we 
have developed generates the transformation very effi- 
ciently. Different numerical techniques applied to com- 
puting the map from arbitrary closed domains to the unit 
disk can be found in Ref.|30|. Our algorithm can be 
shown to be related to the Jacobi method used in these 
studies. 



A. Description of the algorithm 

We define the rough boundary d£ as a single valued 
real function h such that v = h(u) is the equation of the 
boundary. In other words, the boundary is given by the 
parametrised form w — u + ih(u). To comply with the 
framework we have chosen here, we use in the following 
a 27r -periodic h. 

From the particular form of the transformation fi, we 
expand the real and imaginary parts on the boundary d£ 



can be computed by the Fourier transform of hiu). This 
provides a first order approximation of u(x)^ from which 
an improved estimate of Wfe can be obtained, by taking 
the Fourier transform of h(u(x)), i.e. a non-uniform sam- 
pling of the profile h. Iterating this scheme is the basis 
of our algorithm. We will omit for the time being the 
prerequisite on the amplitude of the roughness. We will 
return to this point by considering the stability of the 
algorithm. 

The intermediate quantities appearing at the fcth it- 
eration will be labelled with a superscript (fc). We also 
formulate the algorithm directly in discrete terms suited 
for a numerical implementation. In the remainder of this 
article, all functions will be decomposed over a set of 2n 
discrete values. The number of Fourier modes will thus 
be limited to 2n. We first introduce a series of sampling 
points u^*^'' with j ~ 0, n — I which is initially set to 
an arithmetic series u^^"^ = jT^/n. The sampling of h(u) 
by the u^'^^ gives the array 



(10) 



The discrete Fourier transform of this array is the com- 
plex valued array 



(11) 



m— — n+l 



for — n < j <n. The latter is shortly written as 

aW ^ (12) 

where J- denotes the Fourier transform, which will be 
chosen as the Fast Fourier Transform (FFT) algorithm, 
thus imposing that n is an integer power of 2. The inter- 
mediate mapping lo^^^ is computed from the a^*"'' as 



(fe) 
j 
,(fe) 



r I \ (*-■) 
{i/2n)aQ 



(fe) 



for 



for 



j>0 
J<0 



(13) 



■5Re [Efe^fe.e" 



= 3m [J2k ^ke 



(8) 



If u were equal to x, the second equation would be close 
to a Fourier transform expression of the function h. More 
precisely rewriting the last equation as 



h{u) = 5Re 



. k 



(9) 



we see that the coefficient iUk can be computed from the 
Fourier transform of h{u{x)). The difficulty is that u{x) 
is a priori unknown. However, we note that if the rough- 
ness is small enough, say of order e, u can be written as 
u = X + 0{e). Therefore, identifying x with m is a zeroth 
order approximation. From the latter, the coefficient ujk 



The latter form is obtained from the identification of 
Eq. ^p) and the definition of a^'^^ , taking care of the fact 
that one sum is over positive index, while the other ex- 
tends over the interval [1 — n,n]. Then, one computes 
the series 




Ak) 







L(fe) 



-la 



(k) 



Ak) 



for J > 



for j < 



3 -J -J 3 

This linear transformation is shortly noted as 



(14) 



(15) 



where Q is the above detailed transformation. The form 
of Q is dictated by Eq.JSJi) for positive index, and from 
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the fact that the inverse Fourier transform of b (see be- 
\ow) is real. The new sampHng series is finally obtained 
from 

Jk+i) ^ jn ^ (16) 
■' n 

The equations (| 101 16(1 define one step in the algorithm 
relating w^'^^"'^-' toa;'-'"-'. Briefly we note this step 0^'-'^+"'^^ = 
r(a)('=)). 

The searched function fl is clearly a fixed point of the 
transformation T defined above in a discretized version. 
The uniqueness of the transformation Q results from that 
of the harmonic field in the domain £ with an equipoten- 
tial condition on the boundary and a constant gradient 
perpendicular to the boundary at infinite distance from 
it. Therefore, the only condition to consider is the sta- 
bility of the fixed point. 

B. Stability 

Let us assume that we have an approximate solution 
of the transformation ^{z), from which we compute the 
series Uj . All intermediate quantities computed from the 
exact solution are denoted by a superscript *. Following 
one complete iteration of the algorithm, we obtain the 
following expressions 

Uj — u* + Suj 

hj = h* + h'{u*)5uj 

a = a* +T[h-h*] (17) 

6 = b* + g[a- a*] 

{5u)' = T-^[b-b*] 

where a Taylor expansion of h has been used to estimate 
the values oi h — h* and where indices are omitted when 
unnecessary. The resulting difference {Su)' after one cycle 
is thus 

{5uy = T-^ g T\h:(u*^6uj\ (is) 

Let us introduce the norm 

hii'^Ei^^-i' (19) 

3 

Parseval's theorem relates the above norm in real and 
Fourier spaces according to 

\\T(h)\\^^\\h\\ (20) 

In a similar fashion, the transformation Q does not affect 
the norm 

WQm = m (21) 

Using the two previous results, we can estimate the norm 
of {6u)' as 

< max{\h'\)^\\{5u)f ^ ' 



Therefore, if the absolute value of the slope of the objec- 
tive profile satisfies, 

\h'{u)\ < 1 (23) 

for all u, then the fixed point fl* is attractive for the 
transformation T. It should however be noted that the 
number of modes n should be large enough so that the 
perturbation Su should be small enough to legitimate the 
Taylor expansion of h used in the stability analysis. 

In practice, the convergence is very fast provided the 
sufficient condition H23|l is fulfilled. Moreover it has to 
be noted that one step in the algorithm requires a rather 
limited amount of computing time of order n log(n) (i.e. 
as for a FFT operation). Considering that this compu- 
tation gives the solution of an harmonic problem in a 
semi-infinite domain, this cost appears to be extremely 
low. 

When our algorithm is applied to a simple monochro- 
matic sine (or cosine) profile, h(u) — ^sin(u), it turns 
out that as soon as the condition H23|l is violated (i.e. 
A > 1) the scheme is unstable, and a loop begins to ap- 
pear around the origin where the slope exceeds 1. Thus 
the sufficient condition is also a necessary condition. 

The limit | h' (u) \ < 1 can simply be broken if one uses 
an under-relaxation scheme. The optimum determina- 
tion of the under-relaxation parameter, or the use of 
other algorithms can be found in Ref.j^^ for mapping 
arbitrary domains on the unit disk. The transposition 
of these algorithms to our problem can be worked out in 
details. Other ways to break this limit is to decompose 
the transformation f2 in two (or more) substeps. Sup- 
pose one could map the real axis onto an intermediate 
profile using a first transformation fli and then the in- 
termediate profile onto the objective one using a second 
transformation Q2 ■ The combination of the two transfor- 
mations Q{z) — fl2{^i{z)) is then the searched mapping. 
By breaking the problem into two steps, it is possible that 
each step can be handled by the above presented algo- 
rithm, while the combination of the two giving a profile 
having a slope larger than unity. The difficulty here is 
to device a suited intermediate step. One could consider 
for example filtering the initial objective profile so that 
the filtered profile may fulfill the slope constraint. We 
did not investigate this extension any further. 



C. Convergence and example application 

We present in the following calculations of conformal 
transformations associated with a simple sine interface. 

We will use a norm on the error similar to the one 
introduced in the previous subsection. The distance d 
between the objective profile and the calculated one is 
defined as follows: 

(f ^— \h{u{x)) - n{x)fdx (24) 
Jo 
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FIG. 3: Power spectrum of the function uj computed for a sine 
profile. The three curves corresponds to three amphtudes, 
A =0.3 (•), 0.6 (+) and 0.9 (o). 



FIG. 5: Minimum error obtained using the algorithm de- 
scribed in the text for computing the conformal map on a 
sine profile with amplitude A — 0.5, as a function of the num- 
ber n of modes used in uj. 




FIG. 4; Images of the real axis obtained after the first k 
iterates of the algorithm with the objective profile {h{u) — 
Asin(u)) shown as a bold line. In this particular example 
A = 0.5 and the number of modes is 32. 



It is convenient to make this distance dimensionless, nor- 
malising it by the amplitude of the profile, d* = (d/A) 
where the objective profile has the equation h{u) — 
Asin(u). 

It is worth noting that the problem is far from being 
as simple as it might appear on the surface. In real space 
one single Fourier mode is sufficient to entirely charac- 
terise the interface. The transformation f2 however re- 
quires much more modes. Figure Q shows the power 
spectrum of the u; series, for different amplitudes A. 

The convergence of the algorithm is shown at play on 
Figure where the obtained profile obtained after the 
first few iterations are shown. In this particular example 
A = 0.5 and the number of modes is 32. 

The importance of allowing for enough Fourier modes 
is also illustrated by considering the minimum error d 
obtained as a function of n as shown in Figure 10 . (For 
this particular study we did not resort to a FFT algo- 
rithm to handle any value of n). For A = 0.5 we observe 
that about 20 modes are necessary to reach the single 
precision used in the computation. As the amplitude in- 
creases, the number of modes needed to reach a small 



enough error becomes larger and larger. 



III. SELF-AFFINE BOUNDARIES 

In the description of rough surfaces and interfaces, 
some recent progress has been achieved by recognising 
some scaling invariance properties which has been ob- 
served in a number of real surfaces, and has been shown 
to result naturally in a number of growth models. Recent 
reviewS||5„^, JJ have covered this field. 

Due to the different roles played by the directions nor- 
mal and parallel to the surface, the scaling invariance 
— when applicable — involves different scale factors de- 
pending on orientation, a property called self-affinity. We 
consider here only two dimensional media so that the 
boundary is self-affine if it remains (statistically) invari- 
ant under the transformation 



Ax 
AS 



(25) 



for all values of A. The exponent C is called the "Hurst" 
or roughness exponent. It is characteristic of the scaling 
invariance. From this property, we derive easily that 



{{y{x)^y{x + S)f)=C'd''^ 



(26) 



where C is a prefactor. 

It is noteworthy that the self-affinity property does not 
involve the scaling of any measure. However, studying 
the scaling of the length of the curve, two regimes are 
revealed. For large distances, larger than a scale A, the 
curvilinear length of the profile is simply proportional 
to the projected length along the x axis, hence on can 
identify a trivial fractal dimension equal to 1. On the 
other hand, for distances smaller than A, the arc length 
scales in a non-trivial fashion with the projected length. 
This allows to define a fractal dimension equal to d/ = 
2 — The cross-over scale A between these two regimes 
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is such that the typical slope of the profile is 1, i.e. using 
the notations of Ea. (|26|l . 



(27) 



Once a roughness profile has been measured, a very 
convenient way [sj to check the self-affinity is to com- 
pute the Power Spectral Density (PSD) of the profile. In 
the case of a self-afhne profile of exponent (, the PSD is 
expected to have the following behaviour: 



(28) 



It is important to stress that the approach developed in 
this article is not specific to self-affine boundaries. How- 
ever, being given the practical importance of such bound- 
aries, and the expected generality of scaling results, we 
will essentially focus on self-affine boundaries as practical 
applications of the concepts developed in the framework 
of harmonic field in the vicinity of rough boundaries. 

In view of the form of the transformation Q, and of the 
previous scaling, Ea. H28|l . we introduce a particular set 
of transformation: let us choose 



-1/2-C 



(29) 



where are random Gaussian variable with mean and 
unit variance for the real and imagnary part indepen- 
dently, we can write 



3?e 



dx 



[x + iO) 



= 1 + A. 5m 



/2-Cg-ifca; 



Then for a given set of e^, we can define a maximum 
amplitude such that the mapping is bijective: 



-1 



5m ' ^ 



(31) 



This method gives a short way to generate directly Q, 
transforms which image the real axis to a self-affine in- 
terface shape. This approach is useful to study generic 
properties of self-afhne boundaries. 

When the amplitude A is small enough, u ~ x, and 
thus the series uJk is equal to the Fourier transform of 
the profile. The transformation sends the real axis 
onto a periodic function whose power-spectrum is of the 
form Ea. H28|l . When the amplitude increases, the first 
iteration of the algorithm turns out to be rather approx- 
imative. In order to show that the power spectrum of 
Wfc is not significantly altered by further steps, we show 
in Figure ^ the power spectrum of ujk as compared to 
the initial zeroth order approximation. The result have 
been obtained from an average over 100 profiles with 2048 
modes each. We see that the synthetic generation of the 
transform does not modify the power spectrum of lo co- 
efficients. 

Therefore, we can directly generate mappings lo which 
will transform the real axis into a periodic boundary 
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FIG. 6: Power spectrum of the image of the real axis obtained 
from synthetic transformation Q,, obtained from Ea. l|29|l . The 
roughness exponent used is (" = 0.8. The best power-law 
regression on this power spectrum has a slope s = —1.32, 
which corresponds to an estimated roughness exponent C,' = 
0.82, equal within error-bars to C,. 



which is self-affine with any prescribed roughness expo- 
nent for distances smaller than the period. Such a con- 
struction may appear as artificial in the sense that the 
rough boundary is not imposed but on the contrary it 
results from the choice of the mapping. It is however 
useful, as will be shown later, because it allows to study 
generic properties of harmonic fields close to self-affine 
boundaries. 

The alternative way consists in using the mapping con- 
struction algorithm. We have in a previous section anal- 
ysed the convergence of the algorithm applied to the spe- 
cial case of a sine profile. We now consider the case of a 
self-affine boundary in a similar fashion. This interface 
has been calculated in the real space with 64 modes, and 
we have used 256 modes in the conformal transforma- 
tion. The standard deviation of the height distribution 
is called a. The chosen C exponent chosen for this ex- 
ample is C = 0.8. From Eq.||3U, we note that the max- 
imum amplitude decreases as the number of modes n 
increases. This is natural since as the lower cut-off in the 
scaling regime decreases, the self-affine function will tend 
toward a continuous but non-differentiable curve when 
< C < 1. The distribution of local slopes is indeed ex- 
pected to gets wider and wider as the number of modes 
increases. Quantitatively, Amax ex. n''~^. It is to be noted 
that as n increases, the standard deviation of the height 
a does not increase. It is to be noted that these con- 
clusions are drawn under the hypothesis that the longest 
wavelength remains fixed, here set to 27r. Alternatively, if 
the smaller cut-off and the amplitude of the correspond- 
ing mode were kept constant while increasing the number 
of modes, then the maximum amplitude would remain 
constant. 

As in the previous example (sine profile), we can ob- 
serve on Figure an example of the conformal map 
obtained for a/amax = 0.95, where the maximum stan- 
dard deviation that could be handled by the algorithm 
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FIG. 7: An example of the obtained profile d£ from the con- 
formal map, compared to the objective one, chosen to be a 
self-afBne function with a roughness exponent ^ = 0.8. The 
amplitude of the profile is 95% of the maximum amplitude 
which preserves the convergence of the algorithm. 



without diverging is (Jmax ~ 0.1. We can see on this fig- 
ure that the major differences between the objective and 
calculated profiles occur in areas where the local slope 
is maximum. For roughness amplitude greater than the 
convergence threshold, one can see loops appearing in 
these areas. The convergence speed, the sensitivity to 
the number of modes allowed in the determination of 
the evolution of the minimum error, ... behaved for these 
self-afHne profiles in a similar qualitative way as for the 
simple sine profile. 



IV. GENERIC PROPERTIES OF HARMONIC 
POTENTIALS 

In the following, we show that the knowledge of such a 
conformal transformation allows to solve immediatly har- 
monic problems. We essentially focus here on the case 
where the field is assumed to be uniform far from the 
boundary. This is a typical case as soon as the roughness 
is of small amplitude. This can be seen as an asymp- 
totic expansion focusing here on the small scale details 
of the interface, whereas the matching with the far field 
can be done using a field whose variation is small on 
the scale of the roughness amplitude. We will now focus 
on two problems: perfectly conducting boundary so that 
the potential gradient is normal to the boundary, and 
perfectly insulating boundary where the potential gra- 
dient is parallel to the surface. Since we know how to 
taylor mappings which image the real axis on a generic 
self-afRne boundary, this gives us a key to consider the 
scaling features of harmonic fields in the vicinity of self- 
affine boundaries. 

Harmonic fields are encountered very frequently in na- 
ture. Linear transport involving scalar fields <&, where 
the flux J is proportional to the field gradient plus a 
conservation law in the absence of sources, and in steady 
conditions, div{ J) = 0, imply the harmonic nature of the 



field $, = 0. Heat diffusion obeying Fourier's law, 
gives an harmonic temperature field in steady condition. 
Mass diffusion with Fick's law is a similar example with 
the concentration field. Electric conduction with Ohm's 
law, viscous fiow in confined two dimensional Hclc-Shaw 
cells, vorticity in Stokes flow, ... constitutes a partial list 
of possible applications. 



A. Homegeneous far field 

In this part, for the sake of concreteness, we use the 
case of thermal conduction. We are interested in the 
temperature field T in the region £ limited by the rough 
interface d£. Let us first consider the case of a perfectly 
conducting interface, so that T = Tq for each point of 
the boundary. We impose in the far field a homogeneous 
unit flux of heat. The problem to solve is 

A„T =0 in £ 

T = n on d£ (32) 

VT(w) —fCy if V — oo 

The knowledge of allows then to define 8, image 
field of T in the smooth domain V: Q{z) = T o r2(z) = 
T{w). As A^e = A^T \n'{z)f and fl'{z) ^OinV, the 
resolution of (|32ll in £ is thus equivalent to: 

A^e{z) =0 in V 
y{z) ^To on dV (33) 

VO{z) ^ Ey if y ^ ^oo 

The frontier dV being the real axis, we have imme- 
diately the solution in D: Q{z) = Tq + y and then the 
solution in £ is: 

T{w) =To + ^m[n~\w)] (34) 

Figure (|SJ) shows a set of isotherm curves close to a self- 
affine isotherm boundary. These lines become smoother 
and smoother when the distance to the electrode in- 
creases. The morphology of these isotherm lines have 
some interesting features. If A denotes the distance to 
the boundary, one can observe from the form of the map- 
ping f2 that modes with a wavelength smaller than A will 
be damped whereas longer wavelength modes will only be 
slightly decreased. Therefore, the isotherm curves will be 
similar to the profile up to a low pass filtering. In the case 
of a self-affine boundary, the isotherms will preserve the 
self-affine character with the same exponent, but their 
lower cut-off will increase as the distance to the actual 
boundary, up to the distance of order of the largest wave- 
length. 

Let us now study the temperature gradient. Quite gen- 
erally, we can write the gradient in the complex plane as 

y^Tiw) EE {du + id^a)T{w) = 2d^T (35) 



9 




-0.4 I ' ' ' ' ' ^ 

0.0 2.0 4.0 6.0 

u 



FIG. 8; Example of isotherm curves close to a rough 
self-afBne boundary on which the temperature is constant, 
whereas the temperature gradient is homogeneous and verti- 
cal far from the boundary. 



archetypical problem whose expression is 

A^,T = in £ 
dnT ^0 on dS (37) 
T{w) —fu if f — > — oo 

The solution to this problem can simply be obtained from 
the previous using duality properties of the harmonic 
field. The real part of the previous solution gives the 
answer to the problem. 

T{w) = ne[n-\w)] (38) 

The temperature gradient is then simply 



V^T(zi-)= l + w'(z) (39) 



B. "Equivalent" smooth boundary condition 
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FIG. 9: Profile h{u) of the rough boundary (Symbol •) and 
temperature gradient d„T{u) on the surface of the same 
boundary. We note strong correlations between the two 



From the expression of the temperature field, we have 



(36) 



From the expression of the function to, we see that at a 
large distance from the rough boundary, the term lo' van- 
ishes exponentially. Therefore, one recovers the imposed 
condition for the temperature at infinity i.e. VT — *■ i. 

On Figure we have presented on the same graph the 
profile of a rough electrode and the modulus of the tem- 
perature gradient. One may see quite easily that the field 
is very large (small) in the deepest (highest) areas. This 
field depends naturally both on the local topography and 
on its remote environment. The connection between the 
field and the local topography can be analysed through 
cross-correlations as will be done in a following section. 

The perfectly insulating boundary is the other 



We have seen previously that once we know the confor- 
mal mapping able to transplant the half complex plane 
D onto the rough domain E, we have immcdiatly the so- 
lution of the electrical potential near the rough electrode 
dE. If the roughness amplitude remains below the con- 
vergence threshold, we are now able to solve this prob- 
lem for any kind of boundary. In practice, very often, 
one does not worry about the details of the rough in- 
terface. As we have seen most perturbations die away 
from the boundary exponentially fast. Therefore, know- 
ing the longest wavelength of the boundary gives the scale 
away from the boundary where the field becomes homo- 
geneous. 

This means practically, that if one is interested only 
in the far field, one could replace the rough interface by 
a straight one so that the far field is unperturbed. The 
question we want to address in this section is the follow- 
ing: where should the "equivalent" straight interface be 
located so as to match the asymptotic far field? A zeroth 
order guess is to place it at the geometrical average of the 
height distribution. This will be shown not be an exact 
answer, in the following, we call H the distance between 
the equivalent position and the geometrical average. 

In order to illustrate the problem, let us imagine the 
following experiment. Let us consider an electrolytic 
bath, where the electrical field is homogeneous between 
two opposite electrodes A and B at a, distance L from 
each other. The electrical resistance R of the set-up is 
measured. Then, as illustrated in Figure (|10|) we place 
in the middle of the bath and parallel to the electrodes 
a rough plane C of negligible thickness which is a good 
conductor, so that it can be considered as an equipoten- 
tial. We measure again the electrical resistance of the 
set-up, which is now reduced to i? — AR. What is the 
value of AR? So as to answer this question, we part the 
system in two, A — C and B — C . Each of these two 
problems corresponds to the situation described in the 
introduction of this section. Extrapolating the field from 
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FIG. 10: Schematic illustration of the finite "electrical" thick- 
ness of a rough equipotential. The rough equipotential C 
(shown as a bold curve) is placed in between two remote 
planar electrodes. The presence of the rough equipotential 
reduces the resistance of the medium in a similar way as an 
planar equipotential with a finite thickness. Extrapolating the 
far electrical field from the remote electrodes gives "equiva- 
lent" smooth boundaries shown as dotted lines. Their relative 
distance defines the "electrical" thickness ht. 



electrode A, we find an offset Hi. Similarly from B we 
obtain a different offset iJ2, so that, ignoring the details 
of the perturbed field in the vicinity of C, the rough elec- 
trode will appear as being equivalent to a plane electrode 
of thickness Ht — Hi + i?2- This "electrical thickness" 
has nothing to do with the real thickness of the plane 
considered here to be zero. If the rough electrode has 
the shape of a sine function, of amplitude A and wave- 
length A we will argue below that Ht cx A^ /\. Finally 
it is a simple matter to relate the resistance drop to this 
effective thickness through AR/R — Ht/ L. 

We now revert to the notation of the previous para- 
graph, and deal with the temperature instead of the volt- 
age. For distances away from the rough boundary much 
greater than 27r, (our longest wavelength), all exponential 
terms die out, and hence the far-field can be written 



T(w) « To + '^m(w — LOa) 

where wq is the constant term in the function lo. 
off-set position of the equivalent isotherm is thus 

H = '^m{LU()) 



(40) 
The 

(41) 



Let us first analyse the problem for a small amplitude 
sine boundary of amplitude A and wavelength A. The off- 
set H in the location of the equivalent straight boundary 
is to be normalised by A to obtain a dimensionless quan- 
tity. The latter should be a function of the dimensionless 
ratio A/X. Taylor expansion of this function provides the 
perturbation expansion 



H 



j\2 ^3 

aoA + ai—+ a2TT + OiA^/X^) 
A A^ 



(42) 



A simple argument allows to simplify the latter equation. 
Suppose one would analyse the problem for the profile of 
amplitude —A. The latter is obtained from the former by 



a translation along the x axis by an amount A/2. Thus 
H should be unchanged. This imposes that odd terms in 
the expansion should vanish, hence 



H = ai^+OiAyX^) 
A 



(43) 



Thus the dominant correction is of order A^/X. It can 
be interpreted as the product of the amplitude A and a 
typical slope (A/X). This result holds in the limit of a 
small amplitude and long wavelength. If the wavelength 
goes to zero, clearly the offset should converge to the 
amplitude, but the latter limit cannot be obtained from 
the above Taylor expansion in the small parameter A/X. 

For a sine profile of small amplitude it is possible to 
carry out the computation of the coefficient ai. We 
briefly sketch here the solution. The potential is to be 
computed to second order in A. We revert as in the pre- 
vious sections to a wavelength A — 27r. The solution 
reads 



T{w) = To +3m 



+ 3m 



w — Ae 



— tia _ - —2iw 

2 



(44) 



0{A^) 



The offset can be read from this equation a.s H = 
(1/2) A^. Reincorporating the A dependence, we arrive 
at ai = TT or 



H 



^7r^+ OiA^X^) 
A 



(45) 



This last result is of course only valid for small A values, 
H being bounded by A. On Figures (|lll) and H12|l . we 
can see comparisons between this perturbative calcula- 
tion and the result directly obtained by conformal trans- 
formation. We observe an excellent agreement for small 
amplitudes (resp. large wavelengths) and then the per- 
turbative calculation overestimates H for larger values of 
the amplitude (resp. smaller wavelength). In view of the 
upper bound on the offset and the above perturbation 
expansion, we propose the following form 



2 A . fiT^A 

H « Arctan 

TT V 2A 



(46) 



which fits the data accurately as can be seen on Figs. I|ll|) 
and H12|l . and which reproduces both limiting behaviors 
A ^ and A — > cx). 

The question we ask is now how does the result trans- 
late to a rough profile ? In particular for a self-afhne 
profile, there is no characteristic length scales apart from 
the cut-offs. The product of the amplitude times the 
slope is a scale dependent factor. Is it possible to reach 
quantitative conclusions for such profiles ? 

In order to estimate H for a rough boundary, we use 
the formalism developed for introducing the algorithm. 
We expand the function lo as well as all other interme- 
diate quantities in series of the profile amplitude. Using 
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FIG. 11: Offset H of the equivalent boundary from the geo- 
metrical average of a sine profile of variable amplitude A and 
fixed wavelength X = 2n. The dotted line is the result of the 
perturbation analysis, Eg. 14311 and the symbols are the results 
obtained from a conformal mapping. The dashed curve is the 
proposed fit Ea. ll46ll . 



FIG. 13: Ofi'set H of the equivalent boundary from the geo- 
metrical average of a self-affine profile of variable amplitude 
A, a roughness exponent ^ = 0.8 and 16 modes. The dotted 
line is the result of the perturbation analysis, Eg. 11491 and the 
symbols are the results obtained from a conformal mapping. 



0.00 



The offset can now be expressed as 



— ( J2 hkhk k-J2 hkh k 



\k>0 



k<0 



2n2 



fc>0 



(49) 
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FIG. 12: Offset H of the equivalent boundary from the geo- 
metrical average of a sine profile of amplitude A — 0.05 and 
variable wavelength A. The dotted line is the result of the 
perturbation analysis. Eg. 1431 and the symbols are the re- 
sults obtained from a conformal mapping. The dashed curve 
is the proposed fit Eg. 1461 . 



where hk is the Fourier transform of hj . 

Figure (|13|l gives the evolution of H with the amplitude 
of self-afEne profiles of roughness exponent ( — 0.8, and 
16 modes. Again, we observe that the above given ex- 
pression (|49l) is accurate for small amphtude, but shows 
deviations for larger amphtudes. 

It is interesting to consider the scaling of H observed 
from the generic transformations where u>k are postulated 
to be LUk — Aekk~^^^~'^ . The expectation value of the 
offset {H) reads to dominant order in the amplitude 



the linearity of the transformations Q and J-^ ^\ we 
arrive at 



5fte [oo] 
2n 



= —\"h* 

3 



(47) 



up to third order terms in the amplitude. We now need 
an asymmetric version of Parseval's theorem. Let us com- 
pute the integral for two arbitrary arrays defined in real 
space for u and Fourier space for v 



E"j-^^ '^Mj= 2^EE"j^fee' 

3 3 k 



-ikj 



(48) 



(H) = A- 



-2C 



(50) 



where the extra factor of 2 comes from the expectation 
value of (|ep) = 2, since real and imaginary parts of 
ek are independent Gaussian variables of zero mean and 
unit variance. Depending on the value of the roughness 
exponent C two cases are to be distinguished. 
• For a "persistent" profile — i.e. C > 0.5 — the sum in 
Ea. (|5n|) is dominated by the smallest j, i.e. the longest 
wavelength, and thus the scaling of (H) can be expressed 



as 



(H) 



Z{20- 



27r 



(51) 



where Z{s) = k is the Riemann zeta function. We 
have quitted momentarily the convention that the largest 



12 



wavelength is 27r, hence Xmin and Xmax are respectively 
the smallest and largest cut-off lengths in the profile. In 
this case, A, the amplitude, is such that the largest wave- 
length mode amounts to Aei{27: / \„iax)~^^'^~'' ■ Let us 
introduce the standard deviation of the profile given by 



i/(2n)E^: 



which leads (using Parseval's theorem) to 

^ / 2^ 

T 



(52) 



k 



Z(2C + 1) 



A. 



-2C-1 



(■similar to the scaling Ea. (|26|l . 
Ea. l(5n|) can then be expressed as 



(H) 



Z(2C) 



^(2C + 1) A, 



(53) 



(54) 



The latter equation simply means that the rough pro- 
file behaves as a the simple monochromatic profile. This 
conclusion is however not always valid as is shown in the 
following case. 

• For a "anti-persistent" profile — i.e. C < 0.5 — the sum 
in Ea. f^ is dominated by the largest j, i.e. the shortest 
wavelength, in contrast to the previous persistent case. 

{H)^A'X^<^J^Y~' (55) 

Therefore, we can express the scaling of H in an intrinsic 
fashion as 



A. 



2C-1 



(56) 



In contrast to the persistent case, it appears that the 
offset iJ is dependent on the lower cut-off scale of the 
profile. In fact if Amm is kept fixed, the standard de- 
viation a grows as A^^^. Therefore, one sees that the 
upper scale cut-off disappears, so that H only depends 
on Xmin- In order to see this more clearly, we introduce 
another measure of the roughness which is sensitive to 
the small scale. Let be the norm of the derivative of h: 

= (1/A™a.) / hH^ufdu = (l/2n) V h;[u,f (57) 



which amounts to 



Ar7 



2C-2 



(58) 



From the latter norm, the offset can be written as 

(i/) cx (^2)A„„ (59) 

which is the counterpart of Eq. (|54|) for the antipersistent 
case. 

As a conclusion, the scaling of the offset H is controlled 
by the shortest (resp. longest) scale cut-off of the self- 
affine regime for anti-persistent (resp. persistent) bound- 
aries. 



C. Correlation between local field and topography 

In a preceeding section, we have extracted the expres- 
sion of the temperature gradient as a function of the 
transformation lo. We now use it to investigate the cor- 
relations between the topography and the temperature 
gradient. We study these correlations in the limit of a 
small amplitude. 

The first order perturbation in the temperature gradi- 
ent can be extracted from Ea. H36|) as 



(60) 



where the amplitude of the profile is asssumed to be of 
order e. We introduce the logarithm of the temperature 
gradient denoted ip which can be expressed as 

= log(| VTp) = -2 ^e{J[x)\ + 0[^) (61) 

From now on we will omit the O(e^) term, keeping in 
mind that we focus here only on the dominant term. 

In order to compute the correlation between the gra- 
dient of temperature and the height, we form the cross- 
product and average over u (or x for convenience, since 
their difference is of order e). The expectation value of 
the product is 

((^(w)/i(w)) = -2(3fie[tj'(x)] 3m[w(x)]) = -2^/c|wfc|2 

fc 

(62) 

It is amazing that the same expression appeared when 
computing the offset of the equivalent straight boundary. 

We define now the correlation coefficient a which can 
be identified as the slope of a linear regression between 
h and Its value is 



{ip{u)h{u)) 



- {h\u)) 

since (ft.) = ((/s) = to first order in e. Hence we have 



(63) 



a = -2 



(64) 



This expression holds for any rough boundary of small 
amplitude. 

In the particular case of a self-affine boundary, we as- 
sume as in the previous section that the transformation 
Lo can be taken as the one generated artificially from its 
Fourier decomposition. The latter expression can thus 
be written as 



1-2C 



(65) 



From the latter expression, we have to distinguish be- 
tween persistent ( > 1/2 and anti-persistent C < 1/2 
profiles depending on whether the series is convergent or 
divergent when the number of modes increases. 
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FIG. 14; Average of the logarithm of the temperature gradi- 
ent (if) = (log(|Vrp)), for fixed height h, as a function of h. 
The data points (symbols •) are averages over 1000 profiles 
of small amplitude, with a roughness exponent ( = 0.8, and 
64 Fourier modes. The theoretical prediction is shown as a 
dotted line of slope a*. 



FIG. 15: Evolution of the limit a* as a function of the rough- 
ness exponent (. As approaches 1/2, the coefficient diverges 
as [C, — l/2)~^. The dotted line shows the predicted behav- 
ior and the symbols are the results computed from conformal 
mappings. 



• In the case of a persistent self-afRne boundary, as the 
number of modes increases to infinity, the value of a con- 
verges toward an asymptotic limit a* given by a ratio of 
Riemann Zeta functions 



a 



, Z{2C) 
'Z(1 + 2C) 



(66) 



The divergence of the Zeta function as its argument ap- 
proaches 1, leads to a divergence of a* as {(— l/2)~^. In 
the more general case where Xmax is not set to 27r, the 
above equation should be corrected to 



'Z(1 + 2C) VA 



27r 



(67) 



As a practical illustration of the latter property we 
have studied the correlations between h and ip by aver- 
aging {(fi) at fixed h for 1000 profiles having the same 
characteristics: amplitude A = 0.25Aniax, roughness ex- 
ponent C = 0.8 and 64 Fourier modes. Figure H14() shows 
the evolution of ip versus h. From the Ea. H66|) we es- 
timate a* = -2Z{1.6)/Z{2.6) « -3.50. As shown on 
Figure (|14|l , this value of a* provides an accurate fit to 
the data. The evolution of this coefficient as a function 
of C is shown in Figure l|15|l . 

• The anti-persistent self-affinc profile behaves differently 
from the previous case. The correlation between the 
surface temperature gradient and the height vanishes. 
Mathematically, this result can be traced to the differ- 
ence of behavior of the two series in Ea. (|65|l . However, 
as in the previous section concerning the location of the 
equivalent smooth interface, one can extract the asymp- 
totic behavior of a. 



a 



{<f{u)h{u)) {h'{uf 



= -2 



(h'^iu)) {h{uY 

2 I \min 



1-2C 



27r 



(68) 



This latter result shed some light on the physical mean- 
ing of the previously mentioned divergence. In our pre- 
sentation, we have chosen to fix the largest wavelength 
(set to Xmax — 27r) and amplitude of this mode. In- 
creasing the number of modes implies that the shortest 
wavelength Amm decreases. For roughness exponents in 
the range < C < 1, this implies an algebraic increase of 

with n, while cr^ is bounded. This divergence of i;^ is 
of no importance for the correlation a only if the profile 
is persistent. Otherwise, Ea. (|68|l holds. The perturba- 
tion method used however assumes that both tr^ and 
should be small. The above analysis simply identifies 
which cut-off will dictate its behavior to the correlation. 
The antipersistent case is more suited to the case where 
Xmin is fixed together with its amplitude, while Xmax 
varies. In this case, the a coefficient increases as X]^^" 
as can be read from Eq. (|68|l using the scaling a2 oc A^^^'^^ 
(<r being independent of Xmax)- 



D. Green function for harmonic problems on a 
rough interface 

Up to now, we have only considered harmonic prob- 
lems with a uniform field at infinity. This kind of bound- 
ary condition is of particular interest for problems where 
the scale of variation of the field in the bulk of the solid 
is large compared to the scale of the roughness so that 
an asymptotic development can be performed where the 
matching is to be done on the far field as one focusses on 
the rough boundary. However, from the conformal map- 
ping, one can address more complex types of boundary 
conditions. 

In order to illustrate this, we develop here a particular 
class of solutions which can be used to solve any problem. 
We will consider Green functions which give the field in 
the medium for localized flux / injected in the medium 
from the surface. 
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Let us consider the following problem: a localised flux 
/ = 1 is injected at point (1,0), on the border of the 
unit circle C. The remaining boundary is perfectly insu- 
lating. The same flux is withdrawn at the origin where 
/ = — 1. The harmonic field which fulfills such boundary 
conditions is 



$(z) 



3?e 



log 



(.-1)^ 



(69) 



This potential $ is the Green fmiction for the Domain C. 
Considering the transformation z ^ —i log(2;) maps the 
unit circle to the semi-plane T>. In the transformation, 
the potential $ becomes 



^{z) = - SRe [iz + 2\og (l - e'")] 



(70) 



which is the Green function for a unit flux localised at 
every site {2kiT, 0) for all integer k. At infinity, the po- 
tential approaches ^{z) y. From this Green function 
it is simple to derive the one obtained for a translated 
array of sources. For sources at {xo + 2fc7r, 0), we have 



^z, Xo) = - SRe iz + 2 log (l - e*(^°-^)) 



(71) 



From this latter expression, the Green function for a 
localised and periodic source on the rough profile is ob- 
tained by combining $ and O. The Green function thus 
reads 



(72) 



which gives the potential at point w for a series of sources 
periodically spaced with the same period as the profile 
Wo + 2kn. 



V. CONCLUSION 

We have introduced here conformal mapping tech- 
nique which allows to address harmonic problems in semi- 
infinite domains limited by a rough interface. This map- 
ping is accompanied by an eflicient numerical technique 
which allows us to compute the mapping by a few itera- 
tions of a one dimensional Fourier transform. Moreover, 



this technique provides a natural basis for discussing an- 
alytically some practical applications. 

We then defined and studied the notion of an equiva- 
lent smooth boundary, whose position has been obtained 
exactly in the limit of a small amplitude. This question 
underlined the diff'erences between persistent and anti- 
persistent boundaries, in terms of sensitivity to the lower 
or upper scale cut-off of the self-affine character of the 
boundary. 

We considered the question of correlations between the 
gradient of the harmonic field on the boundary and the 
height of the profile at the same point. The correlation 
has been explicitly computed and shown to converge to 
a precise limit for persistent boundaries. Anti-persistent 
profiles lead to a correlation coefficient which is depen- 
dent on the self-affinity range. 

Finally we have shown that the same mapping can be 
used to address different boundary conditions, including 
the extreme case of a point like source on the boundary 
which is treated exactly. 

Extensions of the above technique are numerous. We 
essentially focused here on static problems involving har- 
monic fields. However, the same mapping may also be 
used in connection with evolution problems such as dif- 
fusion or wave propagation (localisation). Thermal diffu- 
sion in the vicinity of a rough boundary has recently been 
shown to display an anomalous scaling behavior at early 
stages which could be addressed by such methods. The 
a.c. impedence of rough electrodes is another potential 
field of extension which has been studied in recent years. 

It may also be one constitutive brick of a different map- 
ping dealing with different geometries. An example of 
such extensions is the stress intensity factor (i.e. the sin- 
gular behavior of the stress field) at a crack tip. In the 
framework of antiplane elasticity one can compute the lo- 
cal stress intensity factor at the crak tip and relate it to 
the far-field singular behavior. This problem is currently 
being investigated. 

In a companion paper we extend our mapping to 
computations of bi-harmonic fields with applications to 
Stokes flow close to rough boundaries and elastic stress 
fields close to a rough surface. 
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